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Abstract 

In this paper we employ a novel technique combining the Euler Maclaurin formula 
with the saddle point approximation method to obtain the asymptotic behavior (in 
the limit of large representation index J) of generic Wigner matrix elements D'^ M ,(g). 
We use this result to derive asymptotic formulae for the character x J (d) °f an SU{2) 
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1 Introduction 



The saddle point approximation (SPA) is a classical algorithm to determine asymptotic 
behavior of a large class of integrals in some large parameter limit pQ. One uses it when 
exact calculations are either too complex or not very relevant. Recently SPA has been used 
in conjunction with the Euler Maclaurin (EM) formula to derive asymptotic behavior of 
discrete sums [21 E]. In the combined EM SPA scheme corrections to the leading behavior 
come from two sources: the derivative terms in the EM formula and sub leading terms in 
the SPA estimate. 

In this paper we use the EM SPA method to derive the asymptotic behavior of Wigner 
rotation matrix elements. We subsequently use this asymptotic formula to derive the asymp- 
totic behavior of the character of an SU(2) group element. Although our estimate is obtained 
after using twice the EM SPA approximation and once the Stirling approximation for Eu- 
ler 's Gamma functions it turns out to be the exact result. We then proceed to obtain the 
asymptotic expression for Wigner's 3j symbol, recovering with this method the results of [I]. 

Our results are relevant for computing topological (Turaev Viro like [SJ) invariants and 
in connection to the volume conjecture [6j. From a theoretical physics perspective they are 
of consequence for spin foam models [7], Group Field Theory [SJ |9], discretized BF theory 
and lattice gravity [I0],[TT], [12]. Continuous SPA has been extensively used in this context 
to derive asymptotic behaviors of spin foam amplitudes [13], [H], [15], and [IB] . [TTj . [IS] . 

In the recoupling theory of SU(2), the EM SPA method has already been used to obtain 
in a particularly simple way the Ponzano-Regge asymptotic of the 6j symbol [3], [T9]. The 
main strength of this approach is the following. Most relevant quantities in the recoupling 
theory of SU (2) are expressed in Fourier space by discrete sums. In particular, the Wigner 
matrix elements admit a single sum representation [20J. However, generically, the sums 
are alternated hence difficult to handle. Our EM SPA method deals very efficiently with 
alternating signs: generically such signs lead to complex saddle points situated outside the 
initial summation interval. After exchanging the original sums (via the EM formula) for 
integrals one only deforms the integration contour in the complex plane to pass trough the 
saddle points in a completely standard manner. This feature is the crucial strength of our 
method, and allows rapid access to explicit results. The EM SPA method should allow one 
to prove for instance the asymptotic behavior [21] of the 9j symbol. 

The proofs of our three main results (Theorems [TJ [2] and [3]) are straightforward, but 
the shear amount of computations performed renders this a somewhat technical paper. In 
Section [2] we give a quick review of iterated saddle point approximations. In Section [3] we 
establish Theorem [I] and use it in Section [4] to derive the character formula (Theorem [2]). 
Section [5] proves the asymptotic formulae of the 3j symbol (Theorem [3]). Section [6] draws 
the conclusion of our work and discusses the relation between our result for the character 
and the Duistermaat Heckman theorem. The (very detailed) Appendices present explicit 
computations and detail the EM derivative terms. 
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2 Successive saddle point approximations 

We briefly review the iterated SPA approximations. The result of this section justifies the 
use of our asymptote of the Wigner matrices to derive the asymptotic behavior of SU{2) 
characters and Wigner 3j symbols. 

Consider a function / of two real variables. We are interested in evaluating the asymptotic 
behavior of the integral 

dudx e Jf{u ' x \ (1) 

for large J. One can chose to either evaluate I via an SPA in both variables at the same 
time or via two successive SPA, one for each variable. The question is if the two estimates 
coincide. This problem is addressed in full detail in pQ and the answer to the above question 
is yes (for sufficiently smooth functions), with known estimates. Let us give a quick flavor 
of the origin of this result. 

Remark 1. Let f : R X R — > C be a function with and unique critical point (u c , x c ) and non 
degenerate Hessian at (u c , x c ) such that I = f e J ^ u,x ' admits a SPA at large J. Assume that 
the equation d u f(u, x) — admits an unique solution u c = h(x), such that [d^f)(h(x), x) ^ 0. 
Then the SPA of J e J ^ u,x) in both variables (u,x) gives the same estimate as two successive 
SPAs, the first one in u and the second one in x. 



Proof: The simultaneous SPA in u and x yields the estimate 

2tt 



J f(u c ,x c ) 



Jj[dlfdlf-{d u d x ff] 



(2) 



(U C ,X C ) 



The saddle point equation for u, [d u f] (u, x) 
in u yields 




2vr 



dx 



-d 2 j\ 



0, is solved by u c = h(x). Thus a first SPA 

(3) 



,Jf(h(x),x) 



(h(x),x) 



We evaluate eq. ^ by a second SPA, in the x variable. The saddle point equation is 



d 

dx 



f(h(x),x))=[d u f] 



dh 



+ [dxf] 



(h(x),x) 



(4) 



{h(x),x) dx 

and, as [d u f](h(x), x) = 0, the first term above vanishes. The critical point x c is therefore 
solution of [d x f] = 0. The second derivative of f(h(x),x) computes to 

d 2 



(h(x),x) 



dx 

and noting that 
d 



2X f(h(x),x) 



d 

dx V 



{[dxf] 



(h(x),x) 



[d u d x f] 



dh 

(h(x),x) dx 



+ [d 2 x f] 



(h(x),x) 



dx 



[d u f] 



(h(x),x) 



=► [d 2 J] 



dh 

(h(x),x) dx 



+ d x [d u f] 



(h(x),x) 



dh 
dx 



[dm 



(h(x),x) 



(5) 



(6) 



2 



the estimate obtained by two successive SPAs is 



(h(x c ),x c ) 



2tt 

[d u d x f] 2 

[din 



J f(u c ,x c ) 



(7) 



(h(x c ),x c ) 



identical with eq. 

□ 

This remark generalizes PQ, for sufficiently smooth functions of more variables with non 
degenerate critical points. In the sequel we will express the Wigner matrix elements D J MM , 
(up to corrections coming from the EM formula) as integrals which we approximate by a 
first SPA. To compute more involved sums or integrals of products of such matrix elements 
(the character of a SU{2) group element and the 3j symbol) we will substitute the SPA 
approximation for each D J MM , and evaluate the resulting expressions by subsequent SPAs. 



3 Asymptotic formula of a Wigner matrix element 

In this section we prove an asymptotic formula for a Wigner matrix element. Before proceed- 
ing let us mention that many of our results write in terms of angles. In order to avoid issues 
related to the interval of definition of this angles we will always denote them as up = Inw 
for some complex number w with \w\ = 1. 

Our starting point is the classical expression of D J MM , in terms of Euler angles (a, /3, 7) 
in z y z order (see [2"U] ) 

MM ' l,P,7J ~ 4" (J + M-ty.(J-M'-t)\t\(t-M + M')\ 

^2J+M-M'-2t^2t-M+M' ^ 

with £ = cos(/3/2), rj = sin(/3/2). The sum is taken over all t such that all factorials have 
positive argument (hence it has 1 + min{J + M, J — M, J + M', J — M'} terms). We call 
a Wigner matrix generic if its second Euler angle (3 ^ Zti (that is < £ 2 < 1). We define 
the reduced variables x = jj and y = -^7. A priori the asymptotic behavior we derive below 
holds in certain region of the parameters x, y and £ detailed in Appendices [E] and [Cj 

Theorem 1. A generic Wigner matrix element in the spin J representation of a SU(2) 
group element has in the large J limit the asymptotic behavior 

D J xJ , y j(a,(3, 7 ) « e-^-^(-^) l cos[(j + i)0 + xJ^-yJa;-^ , (9) 

with 

A = (i-m 2 -^)-^v^>o, (10) 
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with 0, if) and u the three angles 



1 , 2£ 2 - 1 - xy + 2VA , , ^-x^ + ^VA 
Kp = In — , zip = in 



^(1-^(1-^) v/ea-exi-^ 2 ) 

, ~^ + y£ 2 + VA 
ziu = In — F = (11 

^(i-aa-y 2 ) 

Proof: The proof of Theorem [I] is divided into two steps: first the approximation of eq. (p 
by an integral via the EM formula, and second the evaluation of the latter by an SPA. 
Step 1: In the large J limit the leading behavior of the Wigner matrix element eq. ((8| is 



D J xJ>yJ (a 7 P 7 j) « i- fdu ^K(x : y,u)e J f^ , (12) 



2tt 
where 

f(x, y, u) = —lax — I'jy + vku + (2 + x — y — 2u) In £ + {2u — x + y) In rj 

- x) ln(l -x) + -(l + x) ln(l + x) + -(1 - y) ln(l - y) + -(1 + y) ln(l + y) 
— (1 + a; — it) ln(l + x — u) — (1 — y — u) ln(l — y — u) 

—u \nu — (u — x + y) ln(w — x + y) , (13) 



and 



1 ' y ' J (1 + a; - u)(l - y - u)(u)(u - x + y) ' 1 J 



(15) 



To prove this we rewrite eq. (|8j) in terms of Gamma functions 

^ M/ (a,/3, 7 ) = ^F(J,M,M',t), 
t 

F(J,M,M',t) = e ™t e -™M e ^M> fJ+M-M'^trft-M+M* x 

y/r(j + m + i)r( j -m + i)r( j + M' + i)r( j - m 7 + 1) 
r( j + m - 1 + i)r(j - at - t + i)r(t + i)r(t -m + m' + i) ' 

and use the Euler-Maclaurin formula 

tmax rfmax 

^/i(t) = / /l(t)dt-Si[/l(W)+M*n>ln)] 

f . ^min 

+ E^^^)-^^)], (16) 

where l?2fc are the Bernoulli numbers^] To derive our asymptote we only take into account 
the integral approximation of eq. (15) (the boundary terms are discussed in Appendix [E|, 
hence 

D J MM ,{a,P,y) « J dtF(J,M,M',t) . (17) 



1 Eq. (16 1 holds for all C°° functions hit), such that the sum over k converges. 
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We define u = j hence du = jdt and using the Stirling formula for the Gamma functions 
(see Appendix [A| we get eq. (12). 



Step 2: We now proceed to evaluate the integral (12) by an SPA. Some of the computations 
relevant for this proof are included in Appendix [Bj Denoting the set of saddle points by C, 
the leading asymptotic behavior of a generic Wigner matrix element writes 



D J x j, yJ (a,/3^) 



1 



'27tJ 



;is) 



Our task is to identify C and compute K 



x,y,u* 

dlf)\x,y,u« and f(x,y,u m ). 



The set C. The derivative of / with respect to u is 

d u f = wr — 21n£ + 2lnr] + ln(l + x — u) + ln(l — y — u) — Inn — ln(u — x + y) 
A straightforward computation shows that the saddle points are the solutions of 



(1 + x — u)(l — y — u) 
ou 2 -w[2(l-£ 2 



x 



'1 -£ 2 ) 

' ^ J + u(u - x + y) = 
? 

2 /] + (l-e 2 )(l + x)(l- 2/ ) 



. 



(19) 

(20) 
(21) 



The region of parameters x,y,£ for which the discriminant of eq. (21) is positive gives 
exponentially suppressed matrix elements while the region for which it is zero gives an Airy 
function estimate. Both cases are detailed in Appendix [C| 

In the rest of this proof we treat the region in which the discriminant of eq. (21) is 
negative. We denote by A minus the reduced discriminant, that is 



A 



[i-exe-xy)- 



(x - yf 



> 



and the two saddle points, solutions of eq. (21), write 

x — V rr 



(22) 



(23) 



thus the set of saddle points is C = {«+,«_}. 



Evaluation of f(x,y,u±). We rearrange the terms in eq. (13) to write 
f(x, y, u) = —lax — ijy + (2 + x — y) In £ + (— x + y) In r\ 

+-(1 - x) ln(l -x) + ^(l + x) ln(l + x) + ^(1 - y) ln(l - y) + -(1 + y) ln(l + y) 



(1 + x) ln(l + x — u) — (1 — y) ln(l — y — u) — (—x + y) ln(w — x + y) 
1 -£ 2 (l + x-u)(l - y-u\ 



+u In 



u{u — x + y) 



(24) 



Note that by the saddle point equations the last line in eq. (|24|) is zero for u±. The rest 

(25) 



of eq. (24) computes to (see Appendix B.l for details) 

f(x, y, u±) = —tax — i^y + i[ <fi + xtp — yuo 
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with 



icf) = In 



no = In 



2£ 2 - 1 - xy + 2iVA 
V(l-* 2 )(l-2/T 



= In 



2 



x£ 2 + «V A 



.1^ 



(26) 



Second derivative. The derivative of eq. (19) is 

1 



- 9 u f(x,y,u) 



1 



1 



1 



(27) 



1 + x — u 1 — y — it it -u — x + ?/ 
At the saddle points a straightforward computation shows that the second derivative is (see 



Appendix B.2) 



-dlf) 



1 



m (i-^)(i-y2)e(i-e) 



(aA ± t2y/A~[l + xy - 2£ 2 ] 



(2f 



The prefactor i^. The prefactor K(x, y, u) is 

^/(l- X 2 )(l-y 2 



is: 



w(l + x — w)(l — ?/ — u){u — x + y) 



(29) 



which computes at the saddle points to (see Appendix B.3) 



K 



x,y,u± 



y/(l-x 2 )(l-y 2 )(2^ -l-xy± 2tVA^j 

e(i-a(i-^ 2 ) 2 (i-y 2 ) 2 " 



(30) 



and (30) 



Final evaluation. Before collecting all our previous results we first evaluate, using eq. (28) 

(2£ 2 - 1 - xy ± 2i^[A 



K\ 



x,y,u± 



-dlf)\ 



x,y,u± 



y/(l-x 2 )(l-y 2 )(4A ± z2v / A[l + xy - 2£ 2 ] 

— =- ( 2£ 2 - 1 - ± 2*^ 

v /(l-x 2 )(l- 2 / 2 )(±22 V / A) V 

: (2£ 2 - 1 - xy ± 2*^) 



±z2 v / A v/(l -x 2 )(l -y 2 ) 



(31) 



Comparing eq. (31) with eq. (11) we conclude that 

K\x,y,u± 1 



-«9 2 /)U,^ u± ±t2VA 



(32) 
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Substituting eq. (32) and (25) into eq. (18) we obtain 

-ijax—tj'yy 



D J xJv j(a,f3^) w — ==(—=)"< 

xJ,yJK ,H, U ^J\2^AJ 

^- e t4> e iJ{<t>+ x i>-yu) _|_ \ c -i& e -iJ{<j>+x*i>-yu)\^ ^gg^ 
and a straightforward computation proves Theorem [TJ 



□ 



4 Characters 

In this section we use Theorem [T] to derive an asymptotic formula for the character of a 
577(2) group element. 

Theorem 2. The leading asymptotic behavior of the character of a SU(2) group element 
(with Euler angles (a, (3, 7)) in the J representation, x J (a,/3,7) is 



sin 



(J +1)9 



X>,/3, 7 )^ l — e (34) 

sin 



2 



with 9 defined by 



O P (<* + 7) f o^ 
cos - = cos — cos . (35) 

2 2 2 v ; 



Let us emphasize that up to this point we already performed three different approxima- 
tions: first the EM approximation, second the Stirling approximation and third the SPA 
approximation. To prove Theorem [2] we will use a second EM approximation and a second 



SPA approximation. However, formula (35) is exactly the classical relation between the Eu- 
ler angle parameterization and the 9,n parameterization of an SU(2) group element, thus 
the leading behavior we find (after five levels of approximation) is in fact the exact formula 
of the character! We will discuss this rather surprising result in Section [6j 

Proof of Theorem [2| To establish Theorem [2] we follow again the EM SPA recipe. The 
character \ J of a group element writes 

j 1 
X>,/3, 7 ) = E D mm(«>/3,7) = E D tj,*J^^) > ( 36 ) 

M=-J x=-l 

with x = y the rescaled variable. Note that the step in the second sum is dx = j. The lead- 
ing EM approximation (see end of Appendix [E]) for the character is therefore the continuous 
integral (dropping henceforth the argument (a, 7)) 

X J « J [ dx D J M . (37) 
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as 



We now use Theorem [T] (more precisely eq. (33)) and write a diagonal Wigner matrix element 

(38) 



u xJ,xJ ~ 



4vrJ v / A J 




D Jf(x,x,u + ) 



+ 



-%4> 



,Jf(x,x,u_) 



Note that for diagonal matrix elements the exponents simplify to 

f(x,x,u±) = — i{a + j)x ± 2^0 + x{ip — u) J 



while the discriminant A and angles 0, ip and u> from eq. (11) become 

2£ 2 - 1 - x 2 + 2iv^A 



Kp = In 
iuj = In 



(1 -x 2 ) 
-x(l - £ 2 ) + tVA 



tip = In 



x(l-C)+iVA 

^2 (1 _ a(1 _ = 



ve(i-e)(i- 



A = (i-a(r-^ 2 ) 



(39) 

(40) 
(41) 



We follow the same steps as in the proof of Theorem [T] 
Critical set C x . The derivatives of the exponents for each of the two terms in eq. (38) are 



d x f(x,x,u±) = — i(a + 7) ± — to) ± id x (f) ± ixd x (ip — lo) . 
The derivative of computes to 



\Yi{^e-x 2 + i^l-e? - ln(l - x 2 



2x^\^e 



(1 -x 2 )vF^ 2 



(42) 



(43) 



The difference ip — u computes to 



x(l-£2) +l VA , (Ve-^-^v/l-e 
z(-0 — lo) = In ; ; ^= = m 



-£ 2 ) +zVA 



£ 2 (1 -X 1 



(44) 



and its derivative 



-2x 



^/e -x 2 -txy^i^e 1 



Combining eq. (43) and (45) we have 

d x cj) + xd x (ijj -lo) = 



and the saddle point equations (42) simplify to 



ip — lo = ±(a + 7) . 
Dividing by 2 and exponentiating we get 



(45) 

(46) 
(47) 



Hence the saddle points are solutions of the quadratic equation 



x (1 — 4 J = (4 — x J tan — - — — " 



x 



e 2 sin 2 ^ 



1 - £ 2 cos 2 "pt ' 



Defining a new variable # via the relation cos § = £ cos ^p- the saddle points rewrite 



x 



2 S 2 

2 _£ 2 sin 2 ^ 



sin 



2 6» 



(49) 



(50) 



Taking into account eq. (48) one identifies an unique saddle point (xi) for f(x,x,u + ) and 



an unique saddle point (x 2 ) for /(x,x,u_) with x\ and x 2 given by 

2 



£sin^ £sin^ 



sm 



x 2 



sm 



(51) 



2 2 

Evaluation of the functions and Hessian on C x . Straightforward computations lead to 



2 2 
X l,2 



2N COs2 | 



sm 



2 ' 



A 



asi.a 



2^2 C0S 2 > n 
• 2 9 — ' 

sm 



1,2 



2 2 

Also note that at the saddle points <ft simplifies as 



■ (52) 

sin S 



In 



2£ 2 - 1 - x 2 2 + 2» V /Aj^ 
(l-< 2 ) 



In 



i -asS|-(i-a+¥i-e 2 



In 



2 o @ 6 9 

cos sin — h i2 cos - sin - 

2 2 2 2 



In e 



Substituting the saddle point equations (47) into eq. (39), we see that, at the saddles 
f(xi, xi, u+) = i<p = id , f(x 2 , x 2 , U- 



—u 



-id . 



To evaluate the Hessian at the saddle we first simplify eq. (|42|) using eq. (|46|) hence 
dlf(x,x,u±) = ±zd x (tp — u) = =f2t- 



(1 -x 2 )^ 2 -x 2 



which becomes at the saddle points 

T2i 



cos I 



1 sin 3 f 

T2l i-e cos i 



(53) 
(54) 

(55) 
(56) 



Final evaluation. Using eq. (54) and eq. (56), the SPA of the character eq. (37) is 



1 



X 



2(1 -£ 2 




+ 



sin 3 | 



-10 



-ije 



sin 3 | 



(57) 



which is 



X' 



2 sin 



-i{J+\)e 



sm 



sm 



(58) 

□ 
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5 Asymptotes of 3j symbols 



In this section we employ the asymptotic formula for the Wigner matrices to obtain an 
asymptotic formula for Wigner's 3j symbol. Note that one can use directly the EM SPA 
method to derive this asymptotic starting from the single sum representation of the 3j 
symbol [20]. We take here the alternative route of using the results of Theorem [I] and the 
representation of 3j symbols in terms of Wigner matrices 



dg Di lM ,{g)D J ^ M ,{g)D J * sM ,{g) 



■h 



J\ J2 J3 

Mi M 2 M 3 



Ji J2 J3 
M[ M' 2 M3 



(59) 



where the integral is taken over SU(2) with the normalized Haar measure 



dg 



2tt 



2tt 



da I d'-f dp sin /3 



(60) 



Substituting the asymptote (9) for each matrix element D^i . M ,(g) (i = 1, 2, 3), the integral 



(59) becomes 



dg 



47rJiVA7 



1/2 



47rJ 2 vA^ 



1/2 



1/2 



-iJi(o+7) _ 



Si {^r + Ji (<t>i+Xi'4>i-yiUi)) ^ 



Sri 



(61) 



i=l Si=±l 

We expand (61), perform the integration over a and 7 and change variables from /3 to £ 

• 1 



sin (3d(3 







2 sin — cos —dd = 2 
2 2 H 







d(e) 



n 



to rewrite it as 



E 



=±1 



9 »Ei*(£+/i) 



Sir 



where the index i runs from 1 to 3, ^..7^,0 is a Kronecker symbols and /j is 

ft = Ji [<Pi + x^i - yiUi] . 



(62) 



(63) 



(64) 



We will derive the asymptotic behavior of eq. (63) via an SPA with respect to £ . Note that 
eq. (59) involves two distinct 3j symbols. If one attempts to first set M[ = Mi, and obtain 
a representation of the square of a single 3j symbol, one encounters a very serious technical 
problem. We will see in the sequel that there are two saddle points £± contributing to the 
asymptotic behavior of eq. (63). If one starts by setting = M-, one of the two saddle 
points = 1, and the second derivative in diverges. The contribution of this saddle 
point does not evaluate by a simple Gaussian integration. 

The SPA evaluation of the general case, eq. (63), is a very lengthy computation. We will 
perform it using the classical angular momentum vectors. For large representation index Jj , 
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there exists a classical angular momentum vector Jj in IR 3 of length | Jj| = Jj and projection 
on the Oz axis (of unit vector n) n ■ Ji = Mi. A 3j symbol is then associated to three 
vectors, J\, -h, J3 with |Jj| = Jj and n ■ Ji = Mi = XiJi. By the selection rules the quantum 
numbers Ji respect the triangle inequalities, and M\ + M2 + M3 = 0. This translate into the 
condition that the vectors J, form a triangle Ji + J 2 + J3 = (and n ■ [ J\ + J 2 + J3] = 0). 
The asymptotic behavior of the 3j symbol writes in terms of the angular momentum vectors 
as 

Theorem 3. For large representation indices Ji the 3j symbol has the asymptotic behavior 

U m, Mt ) - T ^-c«[5: W +5)»i + («-/.)*j? + (*- w+g .«») 

y TT{n ■ b) 



with S = Ji A J 2 = J 2 A J3 = J?j A J\, twice the area of the triangle {Ji} and and \I/i. 3 and 
fi ve angles defined as 

^l = ln-' ( ^ A ^ )+a( ^'^ 



S^(nAJ^ 



A3 _ (n A J i) ■ (« A ^3) + in ■ ( J 3 A J, ; ) 



z^ = ln^ v v " y ^ ^, z = l,2. (66) 

^/(«AJ-) 2 (nAJ 3 ) 2 



Before proceeding with the proof of Theorem |3| note that our starting equation ( 59 ) 
involves two distinct 3j symbols. They are each associated to a triple of vectors, Ji, J2, J3 
(| Jj| = Ji and n ■ Ji = XiJ-i) and J[, J' 2) J3 ( | J[\ = Ji, n ■ J[ = yiJi). Note that \Ji\ = \J-\ 
hence the two triangles {Ji} and {J-} are congruent. Consequently there exists a rotation 
which overlaps them. Under this rotation the normal vector n turns into the unit vector k. 
All the geometrical information can therefore be encoded into an unique triple of vectors, 
henceforth denoted J,-, and the two unit vectors ft and k such that \JA = J,-, n ■ Ji = XiJi 



and k ■ Ji = y^i (see figure hj). 




Figure 1: Angular momentum vectors 

Proof of Theorem [3} The proof follows the by now familiar routine of an SPA. We perform 
this evaluation at fixed angular momenta, that is at fixed set of vectors Ji,n, k. 
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The dominant saddle points: The saddle points governing the asymptotic behavior of 

(67) 



eq. (63) are solutions of the equation 



A straightforward computation (see Appendix |D.1[ ) yields 

hence the saddle point equation writes 

= s 1 J l J~K l + s 2 J 2 J~K~ 2 + ssJs^/As . 



(68) 



(69) 



Introducing the angular momentum vectors the saddle point equation becomes after a short 
computation (see Appendix D.2) 

4£ 4 S 2 - 4£ 2 \s 2 + (n ■ k)S 2 - (n ■ S)(k ■ 5)} 

(70) 



+{ [1 + (n • k)] 2 S 2 - 2(n • S) (k ■ S) [l + (n • k)] } = 



for all choices of signs si, s 2 and S3. Dividing by AS 2 , eq. (70) factors as 



.2 l+Jgvfc) 



(L , + (n • S)(k- S) 

1 2 



with roots 



e 



2 1 + (n • fc) 



l + (n-Jfe) (n-S){k-S) 



(71) 



(72) 



2 ' 2 S 2 

again independent of the signs s\,s 2 and S3. To identify the terms contributing to the 
asymptotic of eq. (63) for fixed Ji, ft and k one needs to evaluate J^a/A^ for each of the two 
roots and £ 2 . Using Appendix D.3, we have 

1 



J Z A + 



Ji ■ (n A fc) 



J/Ar 



(5An)(fc-5) + (5A^)(n-S) 



(73) 



To any semiclassical state Ji, n, k we associate six signs, ef and e i defined by 



et-Ji-(nAk) 



I Ji 



(S An)(k-S) + {S Ak){n- S) 



S 2 



Substituting J^a/A^ into the saddle point eq. (|69|) the latter becomes 

1 

2< 



(74) 



(75) 



— * — * — # 

with A + = (n A fc) and A~ 



{SAn)(k-S)+{SAk)(n-S) 



S 2 



As, on the other hand, ^ J; = 0, we 



conclude that at fixed semiclassical state we have two saddle points and two saddle points 
£ 2 contributing 
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The £ 2 saddle point in the term Sj = ef and in the term Sj = — e\ 
The £ 2 saddle point in the term Sj = e~ and in the term Sj = —e\ 



The SPA evaluation of eq. (63) is the sum of this four contributions. 

The second derivative: The derivative of eq. (68) with respect to £ 2 yields 



■(2e - 1 - x lVl 



1 _ £2) SiJi 



2JK; 



(76) 



and the term in the first line cancels (due to the saddle point equation) when evaluating the 
second derivative at the critical points. After Gaussian integration of the dominant saddle 
point contributions, the prefactor in the SPA approximation of eq. (J63|) writes 



(77) 



K 



K — 32 7c 2 s 1 s 2 s 3 i 3 J 1 J2J3\ / A l A 2 A 3 (-df^ ^ s^fi 



The reminder of this paragraphs is devoted to the evaluation of the K for the two roots 
and £ 2 . Substituting the second derivative yields 



K 



± 



-167r 2 sis 2 s 3 z 4 



J1J2J3 



VAiA 2 A; 



(2£ 2 - 1 - x iVi ) 



Taking into account s 2 s 2 s 3 = ef ef, K ± writes 

efefj 2 ^A~fj 3 ^Af [(2£ 2 - l)Jf - J?J?] + Oi 23 



(78) 



-(167T 2 ) 



(79) 



where O123 denotes circular permutations on the indices 1, 2 and 3. Using eq. (72), the 
denominator evaluates, for the £i root, 



1 - (n • kf 



(80) 



while the numerator computes to (see Appendix D.4 for detailed computations and notations) 

1 



e 2 6 3 ^2y A 2 J3W A 3 



( 2 £ l-l)Jl-jnj^ 



+ 



123" 



hence 

Evaluating the denominator in eq. ( 79 ) for £ 2 yields 



(81) 



(82) 



S 2 



l + {n-k) S n S k \fl-(n-k) S H S k 



S 2 
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if/, ^ A2 A S™S k (S H S k ) 2 



S3) 



while a lengthy computation (see Appendix D.4) shows that the numerator is 



e 2 e 3 J 2 y A 2 J31 



(2£ - 1)J* - jfj* 
-S n S k [l - (n ■ kf + 4(n • fc)^- - 4 



+ Ol23 = 
(g"g g ) 2 



J4) 



proving that 



K~ 



-lSn 2 S n S k 



J5) 



Contribution of each saddle: To evaluate the contribution of each saddle point to the 
asymptote of eq. ( 63 ) we first evaluate 



Si [ y + M = Yl Si [{ Ji + \) + x i J i( l ^f) ~ ViJii^t) 



(86) 



Recall that for a fixed semiclassical state only the terms with Sj equal to e^~, —ef, e~ and 
— contribute. We substitute x 3 J 3 = —X2J2 — %iJi an d y 3 J 3 = —y\J\ — V2-h into eq. (86) 
to bring it into the form 



± 



57) 



where 0f , and cjj are the angles 0j, ^ and cjj evaluated at and £ 2 . For each choice 
+ or — in the accolades, one must count both choices of the overall sign. The angles (J)f, 
efipf — efipf, etc. are evaluated by a rather involved computation in Appendix D.5 The 
end results are synthesized below 



le 



ln 



ft ■ ( Ji A S) + iJiS f 



lefipf - zefipf 



k 



ja _ ln (n A Jj) ■ (n A J 3 ) + in ■ ( J 3 A Jj) 



ft A Jj) 2 (n A J 3 ) ; 



J = 1,2 



(88) 



Substituting eq. (88) into eq. (87) yields 



± { E( J * + \) («*« ^ »*B + • + & ■ 

i 

T(k ■ .h)i^f =F (k ■ J 2 )ztt| 3 = ±(Sln T %) , 



2:5 



(89) 
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where denotes 



23 



(90) 



Final evaluation: We put together eq. (82), (85) and (89) and, noting that the two 



contributions form the saddle are complex conjugate to one another we obtain 

Jl J2 J3 \ ( Jl J2 J'3 1 

M 1 M 2 M 3 J \ M[ M' 2 M3 J 

1 1 f P <n H -Q n ) _|_ e -i(n s -%) _|_ ?e i(n s +%) _ jg-iCflfj+ngA _ (g-g 

Taking into account 



Theorem |3] follows. 



(92) 



□ 



6 Conclusion 

Using the EM SPA method we have determined the asymptotic behaviors at large spin J 
of Wigner matrix elements, Wigner 3j symbols and the character x J (fiO °f an SU(2) group 
element g. 

By far the most surprising fact about this computation is that our formula for the charac- 
ter x J (9) is exact. SPA reproducing the exact result for integrals are usually the consequence 
of a Duistermaat Heckman [221 [23] localization property (one of the most famous example of 
this being the Harish Chandra Itzykson Zuber integral [21]). Recall that the Duistermaat- 
Heckman theorem states that a phase space integral 

J n e - lHM , (93) 

where Q is the Liouville form, equals its leading order SPA estimation if the flow of the 
Hamiltonian vector field X = dH) is U(l). To our knowledge all integrals exhibiting 
a localization property (i.e. equaling their leading order SPA approximation) fall in (some 
generalization of) this case. Note that the character of an SU(2) group element can be 
expressed directly double integral by 

x J ( g ) = y^ e HJ,M,t) ~—[ dudx ^/K(x,x,u)e Jf{x > x ' u) + EM. + S. , (94) 

— ' 27T / 

M,t J 

where E.M. denotes corrections coming from the Euler Maclaurin approximation, and S the 
corrections coming from sub leading terms in the Stirling approximation. The double integral 
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in equation ( 94 ) is of the correct form, with symplectic form £1 = ' K(x,x,u)dx A du and 
Hamiltonian f(x,x,u) generating the Hamiltonian flow 



du j u 2 (l + x-u)(l-x-u) f (a+7) (l + x)(l-x-u) \ 

dp ~ V 1-x 2 \ (l-x)(l+x-u)j 

dx lu 2 (l + x — u)(l — x — u) J l7T (1 — £ 2 ) (1 — x — u)(l + x — u) 

dp = V 1 -X 2 ~lf 2 U 2 

Our result can be explained if first, the above flow is U{1) (thus the SPA of the double 
integral is exact) and second the EM and Stirling correction terms cancel, E.M. + S. = 0. The 
alternative, namely that the flow is not U(l) would require an even more subtle cancellation 
of the sub leading correction terms. Either way, the exact result for the character we derive 
in this paper deserves further investigation. 
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Appendix 

In these appendices we detail various technical points and computations. 



A The Stirling approximation 

We detail here the passage from eq. ( Jl7| ) to eq. fll2| . Our starting point is 

D]U,(a,/3,7) ~ J dtF(J,M,M',t) , (A.l) 

with 

F(J,M,M',t) = e^e-^M^M'^j+M-M'-^t-M+M' x 



s/T(j + m + i)r(j - m + i)r(j + W + i)r( j - w + 1) 



(A.2) 



r( J + M - t + l)r(J -M'-t + l)T(t + l)T(t -M + M' + l) 
We use the Stirling formula 

T(n + 1) = n\ w v 7 ^ (-J = V2^e nlnn ~ n , (A.3) 

for all T functions and rescaled variables M = xJ, M' = yj, t = uJ . Collecting all prefactors 
yields 

^ggggl + x)(l- x)(l +y)(l- y) \ a J_ r^f— y (A 4) 

{2nyj\l+x-u){l-y-u){u){u-x + y)) 2nJ V 1 ' V " ' ' K ' 
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and K(x,y,u) takes the form in eq. (14). The "-n" terms in the Stirling approximation add 
to 

\{ ~ J(l + x) - J(l-x) - J(l + y) - J(l-y)} 

- J(l + x - u) - J(l - y - u) - Ju - J(u - x + y) j = , (A. 5) 

which also implies that the coefficient of In J in the exponent cancels. The contribution of 



the T functions eq. (A. 2) is therefore 



|(1 + x) ln(l + x) + (l-x) ln(l - x) + (1 + y) ln(l + y) + (1 - y) ln(l - y)| 



J 
2 

— J|(l + x — u) ln(l + x — u) + (1 — y — u) ln(l — y — u) 
+u ln(-u) + (u — x + y) \n(u — x + y) > . 



The substitution of eq. (A. 6) into eq. (A. 2) yields 



F(J,xJ,yJ,uJ) » ^Ly/K(x,y,u)e J K x *rt 



where f(x,y,u) takes the form in eq. (13), and 



DjU,(a,0, 7 )« / dtF(J,M,M',t)n j dt^j V /K(x,y,u)e Jf ^ 



which reproduces eq. (12) after changing the integration variable to u 



(A.6) 



(A.7) 



(A.8) 



j ■ 



B Evaluations on the critical set. 

In this appendix we present the various evaluations relevant for the proof of Theorem [TJ We 
start by some preliminary computations. Recall that 



A 



1-eW-xy) 



(x - yf 



> 



(B.9) 



As a preliminary we compute the absolute values of the four complex numbers 



u± — x + y = 1 — £ 



2 

x + y 



l + x-u ± = e + X -^Ti^, l-y-u ± = e-^TiVK, (B.10) 



which are 



\ u ^ = (l-f)(l + x)(l-y), 
\l + x-u±\ 2 =£ 2 (l + x)(l + y) 



\u ± -x + y\ 2 = (l-e)(l-x)(l + y) , 
\l-y-u ± \ 2 = e(l-x)(l-y) . (B.ll) 



17 



B.l Evaluation of / at the critical points 



To establish eq. (25) and (26), we evaluate eq. (24) at u± 

f(x, y,u±) = —lax — ijy + (2 + x — y) In £ + (— x + y) In rj 

+-(1 - x) ln(l - x) + ~(1 + x) ln(l + x) + ~(1 - y) ln(l - y) + -(1 + y) ln(l + y) 
— (1 + x) ln(l + x — u±) — (1 — y) ln(l — y — u±) — (— x + y) ln(u± — x + y) . (B.12) 
The real part of f(x,y,u±) is 



$lf{x,y,u±) = (2 + x-y)ln£ + (-x + y) In 77 

+-(1 - x) ln(l _ z ) + 1(1 + x ) ln(l + x) + i(l - y) ln(l - y) + -(1 + y) ln(l + y) 
(1 + x) In 1 1 + x — u±\ — (1 — y) In |1 — y — u±| — (— x + y) In |u ± — x + y| , (B.13) 



and substituting the absolute values computed in eq. (B.ll) yields 
$tf(x,y,u±) = (2 + x-y)ln£ + (-x + y) In 77 

+-(1 - x) ln(l - x) + -(1 + x) ln(l + x) + -(1 - y) ln(l - y) + -(1 + y) ln(l + y) 



(1 + x) 



In 



;-x + y) 



£ 2 (l + x)(l + y) 



(i - v) 



In 



£ 2 (l-x)(l-y) 



In 



(l-fXl-^l+y) 



(B.14) 



Recalling that 1 — £ 2 = 7/ 2 we note that the coefficients of both ln£ and ln(l — £ 2 ) cancel. 
Furthermore, a direct inspection shows that the coefficients of all ln(l — x), ln(l + x), m(l— y) 
and ln(l + y) cancel. Hence 



^f{x,y,u ± ) = . 

Therefore f(x, y, u±) is a purely imaginary number 

1 + x — u± 



(B.15) 



f(x,y,u±) 



-tax — 277/ — (1 + x) In 



|1 + x — u±\ 



(1-y) In 



l-y- U ± 

|1 - y - tt±| 



-x + y) In 



u± — x + y 



\u±-x + y\ ' 
which assumes the form 

f(x, y, u±) = —tax — i-yy ± li^ + xtp — yu 

where the three angles </>, ip and u read 

(1 + x — u + ) (1 — y — u 



(B.16) 



(B.17) 



up = — In 
= — In 
ioj = — In 



|l + x- 


u+\ 


1 + x - 


u + 


|l + x- 


u+\ 


(1-2/- 


U+) 



11 - 



-In 
+ In 
+ In 



11- 


■y-u+\ 


(«+ 


-x + y) 




-x + y\ 


U+ 


— x + y 



\u+ -x + y\ 



(B.18) 
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As the two roots u + and w_ are complex conjugate, one can absorb the various signs in eq. 

(1 + x — U-)(u + — x + y) 



(B.18) and write 



(l + x-u_)(l-y-u_) 
z0 = In , tip — m 



«w = In 



+ £ — — y — M_ 

(1-y -u-)(u + - x + y) 



\l + x — u-\\u + — x + y\ 



|1 — y — U-\\u + — x + y\ 
One by one 0, ^ and to compute by substituting eq. ( B.10[ ) and eq. (B.ll) to 



(B.19) 



i(p = In 



In 
In 



y / e{1 _ x 2 ){1 _ y 2 ) 

e - ^f- + 2^ 2 zVA - (1 - £ 2 )(£ 2 - xy) + 



2£ 2 - 1 - xy + 2^^ 



and 



= In 



ie + gg + ^v 7 ^) (i - e 2 - + 
v/e(i-e)(i-^)(i+y) 2 



(B.20) 



In 



In 



e 2 (l - g) + g - < 2 - ^ - (1 - e 2 )(^ 2 ~ sy) + gjg + t(l + y) VA 

v / e(1 _ e)(1 _ x2)(1+y)2 

-x(l + y)£ 2 + ^ + xy + ^ + z(l+y)VA 



v/£ 2 (i-a(i-^)(i+?/) 2 



n v / e(i-a(i-x 2 ) 



(B.21) 



and finally 



iu = In ■ 



v/e(i-a(i-y 2 )(i-^) 2 



In 



In 



e 2 (i - g) - g + ^ 2 + ^ - (i - m 2 - *y) + gjg: + <l - ^v 7 ^ 

v / e 2 (1 _ a(1 _ y2)(1 _ x) 2 

y(l - x)£ 2 -^ + xy + ^ + z(l - re) VA 



In 



v/e 2 (i-a(i-i/ 2 )(i-^ 



v/e(i-a(i-y 2 ) 

B.2 Evaluation of the second derivative. 



(B.22) 



From eq. (27) we have 

- dlf(x,y,u) 



1 1 

+ 



1 1 

+ 



1 + x — u 1 — y — it u u — x + y 
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(B.23) 



Each term evaluates at the critical points as 





1 


1 + x — u T 


e + ^ ± i 




1 + 


x — u± 


1 + X — U± 2 


£ 2 (l + x)(l- 






1 


l-y + u T 


£ 2 - ± % 


/ — : — 

n/A 


1 - 


y + u± 


\l-y + u T \ 2 


e 2 (i-x)(i- 


-1/) 




1 


u T — x + y 


a-a-- 


=^VA 


u± 


— x + y 


\u T — x + y 2 


(i-e)a- 


x)(l+y) 


1 


u T 








u± 




(l_^)(l + x)(l-y) ' 





The real part of (B.23) is therefore 

e 2 + ^ 



+ 



pi x+y 

S 2 



e 2 (l + .x)(l + y) £ 2 (l-a:)(l-y) 



+ 



(l-e 2 )(l-x)(l + y) (l-e)(l + x)(l-y) ' 



and computes further to 



4A 



(i-x 2 )(i-y 2 )e(i-e)' 



The imaginary part of eq. (B.23) is 

1 



±zVA 



+ 



1 



?(l + x)(l+y) e{l-x){l-y) 
1 1 



-vV 



which finally computes to 



(i-e)(i-x)(i+ y ) (i-e)(i + x)(i- y ) 

1 - 2£ 2 - 



(i-x 2 )(i-y 2 )e(i-e) ' 



B.3 Evaluation of A' 

The prefactor i^U-,y,u ± is 

K 



v/(i-* 2 )(i-y 2 ) 



(1 + x - u±)(l - y - m±)(«±)(m± - a: + y) 



which is, using eq. (B.24), 
K = 



v/(l-a; 2 )(l-y 2 ) 



e(1 _ e)2(1 _ x2)2(1 _ y2) 2 

(i-e 2 



x - y 



and a straightforward computation proves eq. (30). 
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C Real saddle points 

In this section we present the SPA evaluation of a matrix element with 



A = (i-£ 3 )(£ 2 -*v) 



(x - yf 



< . 



(C.31) 



For convenience we denote A' = — A > 0. In this range of parameters the two saddle points 



u± = h±(x;y) = (1 - £ 2 ) + 



^ ' — ± V A' 



(C.32) 



are real. For simplicity suppose that < x < y < 1. A straightforward computation shows 
that <?i_< u + < 1 — hence both roots are in the integration interval. Using the results 
of Appendix B.l the function evaluates at the two saddle points as 

f\ u± = -tax - t/3y ± ($ + x^ - yfi) , 

with 

(2£ 2 - 1 - xy + 2 V / A 7 ) 



$ = In 
= In 
Q = ln 





x 2 )(i 


-y 2 ) 


(-< 2 + 


2 "T" 




veil - 






(£ 2 ?/ - - 


2 -r v 


/AO 



v / e(1 _ e)(1 _ 2/2) - 



(C.33) 

(C.34) 
(C.35) 
(C.36) 



From Appendix |B.2[ we obtain 

-d 2 f 



-4A' =f 2a/A 7 (2^ 2 - 1 - xy) 



«± e(i-e)(i-x 2 )(i-y 2 ) 



(C.37) 



which shows in particular that the maximum of / is w_ (as — <9 2 /| u „ < 0), and the SPA is 
dominated by the latter. In Figure ^ below we represent the function H = $ + x^ff — yVL as 
a function of a; and y, 



The prefactor writes, using Appendix B.3 



V(l-x 2 )(l-i/ 2 )(2e 2 - 1 - - 2y/E)' 



e(i-e)(i-^ 2 ) 2 (i-y 2 ) 2 

hence we get the asymptotic estimate 



(C.38) 



D J x j, yJ (a } P } -f) 



1 



V2%J ^2y/A' 
which indeed is suppressed for large J. 



1/2 



-%a Jx-t-jjy - & -j(®+xy-yn) 



(C.39) 
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Figure 2: The function H = 3> + x\E' — y£l (red) is negative, and vanishes (plane 
z = 0, light blue) when A' (dark blue) vanishes, for £ = 0.1,0.5 and 0.9, from 
the left to the right. 



The case A' = is special. A straightforward calculation shows that under this circum- 
stances 



$ = # = n = o 



(C.40) 



Also, eq. (C.37) implies d^f\ u± = 0. One needs to push the Taylor development around the 
root 



„9 x — y 



to third order 



1 



f(u,x,y) = f +-(u-u o y[d i J] +0(u* 

uo «0 



(C.41) 



(C.42) 



and the Wigner matrix elements has an asymptotic behavior (see [T] 



du v / K(u,x,y)e Jf « e J/|u ° j Ai(a(x, y)[iJ\^)[iJ]~^ + M'(a{x, y)[iJ]^)[iJ]~^\ , (C.43) 

where a(x,y) is some non vanishing smooth real function (determined by K and / evaluated 
at uo, see [T]), Ai is the Airy function of the first kind and Ai' its derivative. At large 
argument the Airy functions behave like 



Ai(C) 



e s 



3 

(2 



The term Ai' is therefore subleading and we have 

du\J K(u, x, y)e J f ~ - 



-Ai'(C) . 



(C.44) 



iJ (ax+-yy- 1 (a(x,y)) § ) 



iJ(a(x,y))V4 



(C.45) 



D Computations for the 3j symbol 

In this appendix we detail at length the various computations required for the proof of 
Theorem |U 
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D.l The first derivative 

To compute the derivative dp ^ s$ («/*), note that d^Ai = — (2£ 2 — 1 — x^). The partial 
derivative of %6i is then 



id. 



2 + 2%-^^- 



(2£ 2 - 1 - + 2i^Ai 



while the derivative of tipi writes 



(2 - i * 2 -j£ iVi ) (2e - 1 - x iVi - 2%JA>) 



(2£ 2 -l-x^) 2 + 4A, 



(D.46) 



-Xi + 1 



1 Xj+Vi 



l-2^ 2 

-XiP+ty/Ei 2^(1-^) 



(2£ 2 - 1 - x^) + *2x iV /A~ 



2 



x^ 2 - «v/A^ 



2^ 2 (i-a(i-^ 2 ) 

We first evaluate 2x;A; - (2£ 2 - 1 - x^i) (2i±M _ ag 

+ ^) 2 



i-2e 
2e(i -a 



. (D.47) 



2a; 
e 2 



-e 4 +r(i+x^) 

Xj(l "I - Xjt/j) Xj ?/j 

' - -''/.H ' ' .' - Cis,). 



4 

Xj + yi 



2£ - 1 - x^ 



/Xj ~T Di 

I— s — ^ 



(D.48) 



hence eq. (D.47) writes 



(l-x 2 )(^-^) «V25 



x^ - 2x 2 e 2 + 2f - 1 



2VA7£ 2 (l-£ 2 )(l-z 2 

2 



(1 - 2£ 2 ) 

2VA7e(i-e)(i-^!) 2e 2 a-a 



2^ 2 

Noting that u(xi,y~i) = ip(—yi, — Xj) the derivative of writes 



(D.49) 



id \^)U) •i 

The derivative of ^ Si(ifi) is then 
% 2 ) X] s *(^) = z 



»(- 



Xi+Vi 



2VA7e(i-e 



(D.50) 



Xj " 



e 2 ^ 



Aj '2VA7e 2 (i-e) 2VA7e 2 (i-e 2 ) 



i 



-2^(1-^ 



2\ I (a^i+j/i) 



- 2£ 2 x^ 



2x/A7e 2 (l-a 
A, : 



(D.51) 
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D.2 The saddle point equation 

We will use in the sequel the short hand notation A B := A ■ B for all vectors A and B. 
Squaring twice the saddle point eq. (69) we obtain, for all signs Sj, 



J 3 2 A 3 - J 2 A 1 - J 2 A 2 



1 2 



4J 1 2 J 2 2 A 1 A 2 . 



We first translate eq. (D.52) in terms of angular momentum vectors 

j-^ = a - e)eji + ejf4 - . 

Then Jf A 3 — J 2 A 1 — Jf A 2 computes to 



(i - e)e 
i 



i 2 



2 J 2 



(J: 



n+k\2 



n+k\2 



+e 



jri jk 
J 3 J 3 



jit jk 

J 1 J 1 



jn jk 
J 2 J 2 



(J: 



n+k\2 



and using J 3 = —J\ — J 2 , eq. (D.52) becomes 



{2(1 - e)efi ■ h + e{J r i4 + 4 4) - \ Jt k Jt k \ 



X 



2(1 - e)( 2 Ji + 2^^^ - ^{Ji +k ? 

2(1 - e 2 )^V 2 + 2£ 2 jfjf - ^(Jf^) 2 



Collecting all terms on the LHS we get 

4(i-m 4 [^ 2 2 -(/i-^ 2 " 

+4(1 - £ 2 )£ 4 [jVPi + JpfJx - Ji ■ h{Jpl + JPi) 

-(i - e 2 )e 2 U 2 (^ +A? ) 2 + JliJ^f - 2/1 • j 2 Jt k J r 2 +u 



H 4 



fc\2 



4jf J* J| J 2 fc - ( Jf J 2 fc + J*J*) 
J? jfiJ^) 2 + J"jf (jf +? ) 2 - (Jf jf + J*J*) J?+*/; 



Eq. (D.56) rewrites 

4(1 -m 4 

-(i-ae 2 



Ji A J 2 



+ 4(1 - e 2 )C n A (Ji A J 2 ) • A; A ( Ji A J 2 ) 



n + fc) A ( Ji A J 2 ) 



1 2 



(nAk) ■ ( Ji A J 2 ) 



-e 2 



n • [(n + k) A ( Ji A J 2 )] fc • [(n + fc) A ( Ji A J 2 )] 



0. 



(D.52) 



(D.53) 



(D.54) 



(D.55) 



(D.56) 



(D.57) 



Using S = J i A J 2 , twice the oriented area of the triangle {Ji}, the saddle point equation 
becomes 







4(l-f) 2 £ 4 S 2 + 4(l-£ 2 )£ 4 

r -* j2 
-2, ,2 l7 * , A o _ ^ 



(n+k)As 



n A S 
(n A k) ■ S 



k AS 

2 



(D.58) 



S-{nAk) S ■ {k An) 
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and dividing by (1 — £ 2 )£ 2 we obtain 

= 4(1 - £ 2 )£ 2 S 2 + 4£ 2 
that is 



has 




kAS 


+ 


S ■ [n A k) 


2 


(n + k) A S 



= 4£ 4 S 2 - 4£ 2 



S 2 + (n- k)S 2 -S H S k 



-S 2 (n A kf +[S A{nAk)] 2 + S 2 (n + kf - (S H + S k ) 



\2 I rin i ak\2 



The last line in eq. (D.60) computes 



-S 2 + S 2 (n ■ k) 2 + (S H ) 2 + (S k ) 2 - 2{n ■ k)S H S k 
+2S 2 + 2S 2 (n ■ k) - (S H ) 2 - (S k ) 2 - 2S n S % 
= [1 + (n • k)] 2 S 2 - 2[1 + (n • k)]S n S n , 



and eq. ( 70 ) follows. 



D.3 Evaluation of JfAf 
Recall that J 2 Aj is 



j 2 A l = (i-e)ej!+ejf^- : ; u 



jn+k\2 



Evaluated for £i 



l+(n-fc) 
2 ' 



eq. (D.62) becomes 



which further simplifies to 

1 



J 2 At = H (n A fc) 2 J- + 2(n • k)JfJ k - ( ) 2 - (J*) 1 
1 



4 



m\2 f jk\2 



4 



(n A k) 2 J 2 + Jf [(n A k) ■ {k A J*)] - J* [(n A fc) • (n A J;)] 
and combining the last two terms this is 

-|(n A k) 2 J 2 + (n A fc) ■ ( J* A (fc A n)) A J t J 

= -J|(n A A;) 2 J 2 + (n A fc) • • (n A fc)) - (n Ak)Jf J, 



hence for £i we get 



Ji- (n Ak) 



n 2 
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Evaluated in £ = ^±|^ - writes 



JA; 



1 - (n ■ k) t S n S k \ (l + (n-k) S n S k ^ , 

J A 



2 J S\ A 2 S 2 
1 + (re ■ k) _ S n S \ jH jk _\( jn+k\2 

2 ^ 4 W< J 



(D.67) 



Combining all the terms common to the RHS in eq. (D.63) and eq. ( D.67[ ), we get 



Ji ■ (n A fc) 



2 £n£fc 
+ 



4S 2 



1 2 



{n-k)J 2 -jfj k 



- J 2 ^ S 



S 4 



(D.68) 



S(Ji ■ (nAJfc)) 



+ 4S n S 



(n A Ji) • (fc A J) 



n afc\2 



-4J? 



£ 2 



But note that S ■ Ji = 0, hence the first term on the RHS above can be written as a double 
vector product, that is 



ja; 



1 J 




4S 12 1 








4S 2 ^ 


: 


1 


AS 4 



Ji A (S A (nAk)) 



+ 4S H S k 



(re A J) • (A? A J) 



4J 



(S n S* 



S 2 



Ji A {nS k + kS H )] 2 -AJ\ 



(S H S k f 



S 2 

2 ■ ^2/ nn ck\2 



S 2 [J A (nS k + kS H ) } - AJ 2 {S H S 
And, as A 2 B 2 = (A ■ B) 2 + (A A B) 2 , we have 

JiAr = 4^1 \S- [Ji A (nS k + kS H )]] 2 + [Ji(S H S k + S k S H )] 2 - Uf(S H S 



(D.69) 



4S 4 

J r [{nAS)S k + (kAS)S H 
AS 4 

D.4 Second derivative 

Using Jja/A/" from eq. (74) and £i, we have 



k\2 



AS 4 



(D.70) 



(2# - l)J( - J?Ji 



+ O 



123 



J 2 nAfe J 3 nAfe [(re A Ji) • (jfe A Ja)] + J 3 nAfe J? A k [(n A J 2 ) • (jfe A J 2 )] 



_ 1 
~ 4 

+^rt«AJ 3 )-(*A J 3 )]} . 
Substituting in the equation above J 3 = — J\ — J 2 , the RHS writes 

' f J 2 " M ' J 2 ™ M '[(re A Ji) • (£ A Ji)] - J" A * Jf A *[(n A J 1 ) ■ (k A J x )\ 

_jHAkjHAk ^ A j 2 ) . ^ A fj] _ jnAkjnAk^ A . ^ A j^] 

(re A Ji) • (A; A Ji) + (re A Ji) • (jfe A J 2 ) 



(D.71) 



_|_ j-nAfc jriAk 



26 



+(n A J 2 ) • (£ A Ji) + (n A J 2 ) • (A? A J 2 )] } , (D.72) 
canceling the appropriate cross terms, the remaining expression factors as 

--{ (nA J!)- J* A *(nA J 2 ) ■ J* h% (k A J x ) - jf A *(£ A J 2 ) } 

= --{nA ((n A fc) A ( Ji A J 2 )) } • jfcA ((n A fc) A ( ^ A J 2 )) J , (D.73) 



developing the double vector products and taking into account that n- (nAk) = k-{nAk) = 0, 
we conclude 



For the £ 2 root we have 
e 2 e 3 J 2 



(2^-l)J 1 2 -« 



+ Oi 23 =-^"5 fe >Afc) 2 



(D.74) 



(2^-l)J 1 2 -(n- J0(fc- Ji 



jSArigk _|_ jSAkgr, 



jSAficjk _|_ jSAkgri 



+ '- 



+ - 



jSAftgk _|_ jSAkCjn 



jSAftcjk _|_ jSAkgf 



S 2 

jSAngk _|_ jSAkgn 



s 2 

jSAngk _|_ jSAkgf 



s 2 



s 2 



Ol 2 3 — 

n A Ji) • (A? A Ji) -2 J 2 
(n A J 2 ) • (A? A J 2 ) - 2 J 2 
(n A J 3 ) ■ (jfe A J 3 ) - 2 J 3 2 



(D.75) 



S 2 J 



S 2 



We substitute again in the equation above J3 = — Ji — J 2 . The coefficient of ^ computes, 
canceling the appropriate cross terms, 



+ 

x 



jSAn gk _|_ jSAk cjri 
jSAn gfc _|_ jSAk gft 



1 2 



(n A Ji) • (A; A Ji) - jf Ari 5 fc + jf Afc S™ (ra A J 2 ) • (k A J 2 ) 



-i 2 



jSAn gk _|_ jSAk 



(n A Ji) • (A; A J 2 ) + (n A J 2 ) ■ (A; A Ji) 



(D.76) 



while the coefficient of 



2S 6 



IS 



Jl 2 



jSAn gk _|_ jSAk gri 



Jo 



jSAri c<k _|_ jSAk q 



SAk rtri 



+2Jx • J 2 



jSAn _|_ jSAH £<• 



SAk art 



jSAn gk _|_ jSAh c< 



SAk an 



(D.77) 



The RHS of eq. (D.75) becomes 
-1 



4S 4 



J 



SAn ofc 



o S K + J^ AK S n (nAJi) 



J$ AH S k + J^ Ak S ri )(kA Ji) 



SAA; 



jf AS 5* + jf A %™) (jfe A J 2 ) 
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S k S n 

' 2S 6 



Ji ( J* AH S k + J$ Ak S n ) - J 2 ( J™ n S* + j?™s 



tSAH ok i jSAk rm 



(D.78) 



which rewrites, combining the appropriate terms into double vector products as 



-1 

49 1 



n A 



[(£ A n) A (Ji A J 2 )^ + (5Afc)A (Ji A J 2 )5"] } 

•{jfe A (S An) A (Ji A J 2 )^ + (5 A fc) A (Ji A J 2 )5"] } 

9^ 9" r ^ ^ ^ ^^-|2 

I (5 An) A(JiA J 2 )S k + (SAk) A (Ji A J 2 )S H , 



2S 6 



and recalling that J\ A J 2 = S the above equation becomes 



(D.79) 



+ 



-1 

49 1 

S k S n 



2S 6 



|nA [(5 An) A55* + (S Ak) ASS'' 
•{jfeA (-9 An) ASS k + (SAk) AS)S H } 

[(SAn) A 5,9* + (5 A fc) A 55™] 2 . 



(D.80) 



— * — * 

We develop the double vector products in the first line and take into account n-(S An) = k- 
(SAk) = 0. For the second line we use (SAA) 2 = S 2 A 2 -(S-A) 2 and S-(SAn) = S-(SAk) = 
to rewrite the equation as 



_ S n S k 



+ 



AS 4 

S k S H 



| - [n- (SAk) S+ (5An)5* + (SAfc)S s 



| - [jfc- (-9 An) S+ (SAn)S k + (S Ak)S H } 
2S4 [(SAn)S~ k + (SAk)S H ] 2 , 



(DM) 



and noting that the cross term in the first scalar product cancel (again as S • (S A n) 
S • (S A k) — 0), and combining the remaining three terms we get 



s n S k 
~4S* 



■S 2 



S ■ (n A jfc) 



2 5*5 



fc OS r 



45 4 



(S A n)S k + (S A k)£T- 



(D.82) 



Factoring S in the second term and using A 2 B 2 = (A ■ B) 2 + (A A B) 2 this rewrites as 



S n S k 



[nAkf 



[nS k -kS H ] 2 [nS k + kS H ] 2 4(S H S 



k\2 



S 2 



+ 



s 2 



s 4 



(D.83) 



thus we conclude 



e 2 e 3 J 2 y A 2 J 3 y A 3 



S H S k 



[n A k) + 4(n • k) 



^ 2 _-l)J 2 -(n- J,)(k- Ji)|+Oi23= 
■ S"5* 4(5^9*) 2 - 



5 2 



5 4 



^D.84) 
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D.5 Function at the saddle points 



We evaluate the relevant angles at the points by substituting eq. (72) and eq. (74) into 



eq. (26) 



D.5.1 The angles <f>f 

For the angles 4>f the direct substitution yields 

(n A Ji) ■ (k A Ji) + iJi Ji - (n A k) 



In 



(nAJiYJikAJ,) 



(n A Ji) ■ (k A f) - 2 J, 



2 S H S k 
i s 2 



J*. 



(5An)5 fe + (5Afc)5 H 



S 2 



^ ^ = In 

^(nA Ji)^(kA Jly 

Consider first the denominator of i(f)~ multiplied by S 2 , namely 

S 2 {n A Ji) ■ (k A Ji) - 2J 2 S H S k + iJ t J t ■ \{S A n)S k + (S A k)S H 
= [SA (n A JI)] ■ [SA (kAJi)] + [S ■ (n A Ji)} [S ■ [k A Ji)} 



-2J 2 S H S k + iJiJi ■ 



(S A n)S k + (S A k)S H 



n ■ ( Ji A S) + iJiS H k ■ ( Ji A S) + iJiS 



hence 



■it. 



i% = In 



n ■ ( Ji A S) + i JiS H 



S\ (n A Ji) 2 



Note that 



(D.85) 



(D.86) 



(D.87) 



n- {Ji A S) + 1 JiS* k ■ (Ji A S) — iJiS k = [S ■ (n A J\)] [S ■ (k A J t )] + J 2 S H S k 



~\~%JiJi 

+j 2 s H s k + ijJi 



S A kS n - nS* 



S 2 (nAJi) ■ (kAJi) ~ [S A(nAJi)] ■ [S A (k A Ji)] 



SA [SA (k An) 



(D.88) 



and developing the double vector products, taking into account S ■ Ji = 0, we conclude 

ze+cf>i = i¥ H - i$\ . (D.89) 
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D.5.2 The angles efipf - efifrf 

We will denote in this section A A B = A aB Direct substitution of and £ 2 yields 

jf - Jf(n ■k)+tf i -(nAk) 



In 



In 



1 - (n • k) 2 ] (ft A Jif 



jf-jf(n-k) + 2jf^ + 



(§An)S k +{SAk)S f 



S 2 



l-(n-k) 2 + A(n-k)^-A 



s 4 



{n A Jif 



To evaluate lefipf — ie 3 ip 3 we take apart the numerator 

jftA(kAn) _|_ % jnAk jHA(kAn) _ 
[nA(fcAn)] -*A [nA(fcAn)] 



jnA(kAn) _ jHAk 
J 3 iJ 3 



J, 



_|_ / jnAk jnA{kAn) _ jfiA(kAn) jnAk\ 



+ Ji-J 3 nA(Un) +JJ^*j; 



i i) ,fc jnAk 



J 

Taking into account n- (k An) = this writes 

-(jfe A n) 2 Jf J 3 ™ + J x • J 3 (£ A n) 2 + iJx 
= (k A nf(n A J x ) ■ (n A J 3 ) - iJ x 

hence 



J?. A 



(ft A k)A(n A (k An)) 



J 3 A 



n(nAk) 2 \ , 



(nA Ji) 2 (nA J 3 ) 2 
To evaluate ^e^i* — ^3^3 we again take apart the numerator 

$An^t+SAk 



/ nA(fcAn)+2n^f- > 



SAn^+SAk^-^ { jnA(kAn)+2n^§ h 



lJ 3 ' ^" 



The real part is 



rcA(fcAn)+2n^f^ nA(fcAn)+2n^f^ 

J 3 f J x 



SAnfz+SAkfz SAnfz+SAkfz 
^3 



(D.90) 



(D.91) 



(D.92) 



(D.93) 



. (D.94) 



, [nA(fcAn)+2n^f^] -*A [nA(fcAn)+2n^|^] 



+ Jl-J'c 



+ 2„^ 



which rewrites, taking into account S • Ji — 0, 

on cfc > , Cm qk 

nJ* AH — (k A ft) J? + 2 J x A n — — I • ( nJ 3 fcA " — (k A n)J* + 2J 3 An 



S 2 



S 2 
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qk _ qit 



.5* 



+Ji ■ -h 



[ft A (k An)] 2 + 4 



' ^ ( J3 g 2 + ^3 ^2 y 

{S H S k ) 2 



s 4 



On . 
S 2 
+ S 2 



Developing the products in the first line we get 



/ on ok\2 
2 trim a\° ° 



_ jkAn jfcAn _ ^ A j^y jn jn _ 4 



S 2 



s 4 



( Ji An) ■ (J 3 An) 



) S H S k 
' S 2 



( J! A n) • (Jfe A n) Jf + ( J 3 An) ■ {k A n) jf 



^2 



Jf S fc + J*5^) ( J 3 "5 fc + J 3 fc 5" ) + Ji • J 3 ( (n A k) 2 + 



(nS fc + kS H ) 2 
S 2 



which is, expanding the second line, 



(D.96) 



(D.97) 



( Ji A n) • ( J 3 A n) ( (n A k) 

S n S k 



n ok\2 



(S n S 



s 2 



jkAri jkAn 



+2- 



S 2 
+Ji ■ J3- 



Jf J 3 n + J^Jl - 2(n ■ fc) J?J£ 
~S 2 ' 



S 2 



J?S k + J k S H J*S k + J k S H 



(D.98) 



Combining the cross terms in the second line and using J^ hn J^ hn = J™ Ak J 3 fe we obtain 



(Ji A n) ■ ( J 3 An)[(n Ak) 2 - A 



it ok\2 



(S n S 



S 2 



4{n-k)J?J. 



S H S k 
S 2 



1 



jnAk jitAk x T SA(tiAk) T SA(nAk) , -, 
J 1 J 3 - —J l J 3 ,/j • J: 



(nS k + kS H ) 2 
S 2 



(D.99) 



and computing the middle term on the second line taking into account 5* ■ Jj = 0, we obtain 



(JiAn) - (J 3 An) (nAfc) 2 -4 



it ok \2 



(S n S 



S 2 

T r (nS k + kS H ) 2 - (nS k - kS H ) 2 

+ <A " ^3 7^ j 



4{n-k)J?4 



S n S k 
S 2 



S 2 



hence the real part is 



(Ji A n) • (J 3 A n) \(n A k) 2 - 4^^£ + ' k) 



S 2 



S 2 



The imaginary part of the numerator ( |D.94 ) writes 



SAn^j+SAk^ nA(kAn)+2n^- nA(kAn)+2n^- SAnj^+SAk^ 

J l J 3 — J l J 3 



We start by expressing it as 



A 


nA(kAri)+2n^§ k - 


A 


SAit^j+SAkj^ 






■Js 





(D.100) 



(D.101) 



(D.102) 
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Shn^+Shk^ 



(D.103) 



as the cross terms in the development of the two scalar products cancel. This computes 
further to 



»i A n 

+ (nJ kAH - (jfe A n)4 l + 2J 3 A n 
Grouping together similar terms we get 



(nj^n _ {kA n)J n + 2Ji A . 5 ^n_ + jk_^ 



S 2 

S n S k 
S 2 



Qk 

S 2 



(D.104) 



£2 



(j 3 A J x ) • ((£ An) Anj + ^"(^ A Ji) • ((£ A n) A fc) 



■ [(A; A n) • 5"] — ( J 3 A Ji ) ■ (n A k) 



(D.105) 



n\2 ok 



+2^i[(nA(J 3 A Ji)) An] • 5 + 2^^- [(fc A ( J 3 A Jx)) An] -S. 
Recognizing S — J 3 A Ji, the above writes 

qnqk _ fS'"'") 2 - S*™ 

— [5«(fl. fc) - + ^L[5« - S fe (n • fc)] + -[5- (n A A;)] 2 



+2 



S 4 



/ qn\2 qk 

2 { L ( -S H S k + {n-k)S 2 



= S n 
= S n 



(ft A k) 2 — 



+ 



[n A kf + 4 



(nS k - kS n ) 2 (S k ) 2 , (5 



->n\2 



S 2 



-4 



(S*\S fc ) 2 (S k ) 2 n _ A 



^S 2 " J 
(D.106) 



S 2 S 4 
In conclusion le^ipi — le^ip^ is 

(nA Ji)-(nA J 3 )+«n-(J 3 A Ji) 



^ ^ — ?e 3 ^3 = In 

^(nA Jx) 2 (nA J 3 ) 2 

Following similar manipulations we get 

(n A J 2 ) • (n A J 3 ) + m • ( J 3 A J 2 ) 



13 



(D.107) 



tef^f - tefipf = In 
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(nAJ 2 ) 2 (nAJ 3 )' 



(D.108) 



For the angles Ui, recall that can be written in terms of Note that due 

to the choice of the determination of the yAf the correct relation is u>t- = —il> + - w and 

v n,k —k,—n 

uZr = il>~r Moreover, as ^ r = — we conclude 

n,k —k,—n K K 



13 

k 



ie 



(D.109) 
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E Boundary terms in the Euler Maclaurin formula 

Using the short hand notation F(t) for F(J, M, M' ,t), the remainder terms in the EM 
formula write 

- b x [F(t max ) + F(t min )] + E pyy - i r(2fc_1) (*mi n )] • (E.no) 

In this section we deal with generic Wigner matrices, that is we consistently assume that 
< £ 2 < 1. Note that 

t min = max{0, M - M'} , t max = min{ J + M,J- M'} . (E.lll) 

For simplicity we will detail the diagonal matrix elements M = M' . By continuity the region 
in which our results apply extends to some strip \M — M'\ < P. For such elements i m i n = 
and, for M > 0, t max = J — M. The Stirling approximations become easily upper and lower 
bounds, at the price of some constants, thus by Appendix [X] we obtain 

■^/K{x,x,u)e mf{x ' x ' u) < \F{t)\ < -—^/K{x,x,u)e mf{x ' x ' u) , (E.112) 



with 



f(x, x, u) = —i(ct + 7)x + iTcu + (1 — u) ln£ 2 + u ln(l — £ 2 ) 
+ (1 - x) ln(l - x) + (1 + x) ln(l + x)- 2u\nu 

-(1 + x-u) ln(l + x - u) - (1 - x - u) ln(l - x - u) . (E.113) 



and 



K (*> u ) = 7TT { \rr^ V~2 ■ (E.114) 

(1 + x — u)(l — x — u)u A 

The behavior of the higher derivative terms in the EM formula is governed by F^ k Ht m m ) 
and F^ k \t ma _ x ). To see this, collect all factors depending on t in F(t) and write 



F(t) = q{J,M) p JiM {t) 

PJ,M(t) 



T(J + M-t + l)r(J - M - t + l)[T(t + l)] 2 ' 
A:=i[Tr-2\n£-2\nr)}. (E.115) 



Hence F^ = q(J, M)p JM , and the first derivatives expresses in terms of 

d r , r'(j + M - t + 1) 



dt rJ,lvlKJ " J ' Mw r 1 r(j + M-t + i) 

ry-M-t + i) r'(Hi) | 
r(j-M-t + i) r(t + i)J 

= p JM (t){A + ^°\j + M -t + 1) 

+^°\j — M — t + 1) — 2^ 0) (t + 1)} , (E.116) 
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Figure 3: Shaded region where the EM corrections are exponentially suppressed 



with -ip(°\t) denoting the digamma function. For integer arguments 

m 1 

^(o) (m + 1) = _ 7o + ^ 



k=l 



(E.117) 



hence < C \nJ \F(t)\ for some constant C. Higher order derivatives of eq. (E.116) 

write in terms of higher order polygamma functions ip {n) = d n ^ 0) /dt n . For all k, tp {2k) (X) < 
ip(°\X) at large X, therefore the fc'th derivative is dominated by 



2fc-l 



(E.118) 



i=i 



Then < C(ln J) k \F(t)\ for some constant C. 



From eq. (E.112) we conclude that both |F(t min )| and |-F(t max )|, as well as all their 
derivatives are a priori exponentially suppressed in the region where 3ftf(x, y, w m in) < and 
3lf(x, y, w max ) < 0. As 



$tf{x,x,0) 
$lf(x, x, 1 — x) 



ln^ 2 

xln^ 2 + (l-s)ln(l -£ 2 ) 

+ (1 + x) ln(l +x)-(l-x) ln(l 



x 



2xln(2x) , (E.119) 



we conclude that the derivative corrections coming from t m ; n = are always suppressed term 
by term. However the situation is markedly different for the corrections coming from t max = 
J — M. At fixed £ 2 , the corrections are exponentially suppressed for x close enough to either 

or 1, but the maximum of 3lf(x, x,l — x), achieved for x — — — 



is 



_ . ln^p^>0, 

hence there exists some interval in which, term by term, the derivative terms are bounded 
from below by an exponential blow up. In this region our EM SPA approximation should a 
priori fail (see also figure [3]). 



C 



A second set of EM derivative terms come when passing from eq. (|36|) to eq. (|37|) 
involving derivatives jjfycDij tX j\ x =±i- Using Appendix 
derivatives yield some function times D J xJxJ . As 



D" 



j,- j 



is) = e J e + ^ )J , D J jj{g) = e J e-^ J 



eq. (|C.39|) we note that all these 

(E.120) 



all such derivative terms are exponentially suppressed for large J . 
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